Giant Nernst Effect due to Fluctuating Cooper Pairs in Superconductors 
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A theory of the fluctuation-induced Nernst effect is developed for arbitrary magnetic fields and temperatures 
beyond the upper critical field line in a two-dimensional superconductor. First, we derive a simple phenomeno- 
logical formula for the Nernst coefficient, which naturally explains the giant Nernst signal due to fluctuating 
Cooper pairs. The latter is shown to be large even far from the transition and may exceed by orders of mag- 
nitude the Fermi liquid terms. We also present a complete microscopic calculation (which includes quantum 
fluctuations) of the Nernst coefficient and give its asymptotic dependencies in various regions on the phase 
diagram. It is argued that the magnitude and the behavior of the Nernst signal observed experimentally in 
disordered superconducting films can be well-understood on the basis of the superconducting fluctuation theory. 

PACS numbers: 74.40.+k, 74.25.Fy, 72.15.Jf 



A series of recent experimental studies has revealed an 
anomalously strong thermomagnetic signal in the normal state 
of the high-temperature superconductors [ HBIsl 0> El H S 
|5[] and disordered superconducting films 1 1 Qui 1 [I . In the pio- 
neering experiment | U, Xu et al. observed a sizeable Nernst 
effect in the La2- x Sr f Cu04 compounds up to 130 K, well 
above the transition temperature, T c . This and further similar 
experiments on the cuprates have sparked theoretical interest 
in the thermomagnetic phenomena. Theoretical approaches to 
the anomalously large Nernst-Ettingshausen effect currently 
include models based on the proximity to a quantum critical 
point 01211 . vortex motion in the pseudogap phase [2. 
as well as a superconducting fluctuation scenario |L5| 
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While the two former theories are specific to the cuprate su- 
perconductors, the latter scenario should apply to other more 
conventional superconducting systems as well. Very recently, 
a large Nernst coefficient was observed in the normal state of 
disordered superconducting films [Till . These supercon- 
ducting films are likely to be well-described by the usual BCS 
model and, hence, the new experimental measurements pro- 
vide indication that the superconducting fluctuations are likely 
to be the key to understanding the underlying physics of the 
giant thermomagnetic response. 

Various groups have previously calculated the fluctuation- 
induced Nernst coefficient in the vicinity of the classical tran- 
sition lll5[ll6lll7lll8Ul9ll . However, these analyses were lim- 
ited to the case of very weak magnetic fields and temperatures 
close to the zero-field transition, when Landau quantization of 
the fluctuating Cooper pair motion can be neglected. In exper- 
iment, however, other parts of the phase diagram (in particular 
strong fields) are obviously important and how the quantized 
motion of fluctuating pairs would figure into the thermomag- 
netic response has remained unclear. In this Letter we clar- 
ify this physics, explaining the origin of the giant fluctuation 
Nernst-Ettingshausen effect, and develop a complete micro- 
scopic theory of Gaussian superconducting fluctuations at ar- 
bitrary magnetic fields and temperatures. 



We start with a qualitative discussion of the Nernst-Ettings- 
hausen effect. Consider a conductor in the presence of a 
magnetic field, H z , and electric field, E y , directed along the 
z- and i/-axes respectively. The charged carriers subject to 
these crossed fields acquire a drift velocity v x = cEy/H- 
in the x-direction. The latter would result in the appear- 
ance of a transverse current j x = neu x . When the cir- 
cuit is broken, no current flows, and the drift of carriers is 
prevented by the spacial variation of the electric potential: 
V x (p = -E x = {nec/o-)(E y /H z ), where cr is the conductiv- 
ity. Due to electroneutrality, this generates the spacial gradi- 
ent of the chemical potential: V x fi(n, T) + eV x <p = 0, which 
corresponds to the appearance of the transverse temperature 
gradient V X T = V x p(dp/dTy l along the x-direction. Hence, 
the Nernst coefficient can be expressed in terms of the full 
temperature derivative of the chemical potential: 



v N 



dp 



{-V X T)H Z ne 2 c dT 



(1) 



E.g., in a degenerate electron gas, the chemical potential 
p(T) = po - (7T 2 T 2 / 6)(d In v/ dp), where v(p) is the density 
of states, and one easily reproduces the value of the Nernst 
coefficient in a normal metal: v N = (n 2 T /3mc)(dT / 'dp) If20l 
2111 . where t is the elastic scattering time (here and below h = 
ks — 1). Thus the Nernst effect in metals is small due to the 
large value of the Fermi energy. 

The simple form of Eq. ([TJ, which to the best of our knowl- 
edge has not appeared previously in the literature, suggests 
that in order to get a large Nernst signal, a strong tempera- 
ture dependence of the chemical potential of carriers is re- 
quired. This can be achieved in the vicinity of the supercon- 
ducting transition where a new type of carriers (fluctuating 
Cooper pairs) appear besides the normal electrons. These ex- 
citations are unstable, have the characteristic lifetime of order 
tgl = 7r/8(r - T c ), and form an interacting Bose gas with a 
variable number of particles. In two dimensions, their concen- 



tration is n^lfT) = (mT c /n)]n[Tc/(T — T c )] II22H . In the vicin 



ity of T c , the chemical potential of fluctuating Cooper pairs 
can be found by identifying its value in the Bose distribution 
to give «c 2 p.(T) above. This leads to nf^.{T) = T c - T, Since 

(21 

dfi\L JdT = -1, the fluctuation contribution to the Nernst sig- 
nal (fJI exceeds parametrically the Fermi liquid term. In this 
sense it is similar to the fluctuation diamagnetism (which also 
exceeds the Landau/Pauli terms and is effectively a correction 
to the perfect diamagnetism of a superconductor). Based on 
Eq. (Q3 and using the known expression forparaconductivity 
in a magnetic field, cr^ - (e 2 /2e)F(e/2Ji) 022I1 . one finds (with 
logarithmic accuracy) the value of the Nernst coefficient due 
to fluctuating Cooper pairs: 





FIG. 1: The Aslamazov-Larkin (AL) and density-of-states (DOS) di- 
agrams for the thermoelectric response p xlJ . The DOS diagram has a 
symmetric counterpart. The white and black circles correspond to 
the different heat and electric vertices, the shadowed blocks repre- 
sent cooperons, and the wavy lines denote the fluctuation propagator 
(see text). All objects on these graphs are generally matrices in the 
Landau basis. 
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[mc(J-T c )Y\ x»l, 



(2) 



mc T -T c I (meDH) 1 , x «: 1 , 

F(x) = x 2 [^(1/2 + x) - <Hjc) - l/(2x)] , (3) 

where x = e/2h, e = ln(T/T c ) and h = H/H C 2(0) are the re- 
duced temperature and magnetic field, H C 2(0) = 4cT c /neD is 
the linearly extrapolated value of the upper critical field, and 
D is the diffusion coefficient. The estimate © corresponds 
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Aslamazov-Larkin (AL) term, it dominates in the Ginzburg- 
Landau (GL) region in the immediate vicinity of the classi- 
cal transition and competes with the negative density-of-states 
(DOS) contribution everywhere else. It is interesting that the 
hierarchy of the DOS diagrams in the problem under consider- 
ation in unusual: Only the graphs containing three cooperons 
(see Fig. Q]) turn out to be important. The AL and DOS con- 
tributions, and the fluctuation magnetization can be written as 



to the results of Ginzburg-Landau (GL) II15LI18H and diagram- n xy (,i\- A V „T V V a x CO, \i (OA 

. rr-\, ., ... , , , ., .. ^pjy^v) - *VHl / / Hmn^nrnHnm-Dnm L m( il k+v)^n(±*k), 



matic lfl6ll treatment in the vicinity of the classical transition 
To develop a quantitative theory of observable thermoelec- 
tric transport, we have to recall a deep relation between the 
fluctuation Nernst effect and magnetization as emphasized in 
Refs. Il5l l23l 12411 : In the presence of a magnetic field, the 



measurable transport heat current differs from the micro- 
scopic heat current j G by the circular magnetization current 
= cM x E, where M is the induced magnetization. As a 
result, the thermoelectric tensor B al3 relating j® a = TB^E 13 
with the applied electric field E can be found as a sum of the 
kinetic, /3 a/3 , and thermodynamic,/^, contributions: 



M • 



0* 



c afiy 



cM 7 IT. 



(4) 



The term ft" 13 can be expressed via the thermal and quantum 
mechanical averaging of the electric-heat currents correlator 
g^(Wy) = (/ ff ( _w v)7 G ^(w v )) by analytic continuation to real 
frequencies: ]?$ = T~ l lim^o Im G a/J ( —iio + 0)/ll>, while the 
term ffj^ accounts for the magnetization heat current 
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The general expression 122l 12511 for the fluctuation mag- 
netization M(T, H) has been calculated previously in the GL 
region 122112611 and at low temperatures, close to H C 2(0) 12711 . 

Now we proceed with the microscopic calculation of the 
Nernst coefficient v N (T,H) = R D B xy /H at arbitrary T and H. 
We concentrate on the calculation of the linear response oper- 
ator Q xy (a) v ) and its analytic continuation to real frequencies. 
In order to get the most general expression we follow Ref. 12711 
and perform calculations in the Landau basis, without expand- 
ing the Green's functions, propagators, current and heat ver- 
tices in the magnetic field. Within this approach, the fluctu- 
ation part of the correlator Q xy (u> v ) is generally represented 
by ten diagrams 1221 12711 . However, in the case of the Nernst 
effect, the Maki-Thompson contribution to Nernst coefficient 
can be shown to be exactly zero. What concerns the positive 



fit n,m 



(5) 



2<2 r ?o S (*>v) = 4v w r J] J] q^n^qlM^x (6) 



fit n,m 



M^-^VHT^lnZ^). 



(7) 



fit n 



HereL„(Q,) = -v" 1 [ln(T/T c ) + MW) ~ tKl/2)] -1 is the 
fluctuation propagator, if/„{Q) is a short-hand notation for 
i/r[l/2 + (Q + a„)/4nT], with a„ = (4eDH/c)(n + 1 /2) being 
the Landau spectrum, vh = eH/nc, and the matrix elements 
of the momentum operator in the Landau basis are given by 

The three-Green-function blocks with two cooperons and 
electric (—e) or heat (i [e/ + e; +v ] /2) vertices, B„„ (f2^, u> v ) 
and S, ( £' (fib ^V)> and the six-Green-function block with three 
cooperons and electric and heat vertices, S^P(Oj, lj v ), shown 
in Fig. Q] should be calculated for an impure metal exactl y (th e 
important ingredients can be found, e.g., in Refs. 1221 12711 '). 
Performing summation over the fermionic Matsubara energies 
we find the following explicit expressions (oj v > 0): 



B^hQ. k ,0J v ) = evD 



i/f m (o) Y + \Sl k \) - i[f„(fik\) 



oj v + a m - a„ 
i/r n (a) v + lOfc+yl) - ^ m (|£2jk +v |) 



-ivD 



a> v - a m + a„ 



, (8) 



x 



(Ot - a m )t//,„(\D. k \ + cj v ) - (Q-k+v - a„)4i n {\D. k \) 



to v + a m - a„ 
(D. k+V + a n )i/f n (\£l k+Y \ + co v ) - (Cl k + a m )i// m (\Q. k+v \) 



a> v + a n - a„ 



(9) 
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X%®(£l k ,(O v ) = -ievD 2 



cj v - a m + a„ 



— l//„(\Q. k+ y\ + LOy) 



D. k - a„ 



{u v + a m - a n ) 2 
Q. k + a, 



(i[s m (\n k \ + <o v ) - iM|Q*|)) 



(a> v - a m + a„) 



- (i/f n (\o. k+Y \ + a> v ) - iff,„(\n k+v \)) 



(10) 



The general calculation of Eqs. ©-([Til is straightforward 
but very cumbersome. However, one can identify nine quali- 
tatively different regions on the phase diagram (Fig.[2]i, where 
the asymptotic behavior has a simple analytical form. The 
corresponding nine asymptotes are presented below. Before 
proceeding to the details, we point out that the resulting ex- 
pression for the Nernst coefficient appears to be universal: 
Contrary to the behavior of fluctuation conductivity 12211 . the 
function /3 xy (T,H) depends only on T/T c and H/H c2 (0) and 
does not depend on the elastic scattering time t. This is the 
magnetization contribution, fy^, which regularizes the other- 
wise divergent (and thus r-dependent) terms in p xy . 

We start with discussing the classical regime close to the 
critical temperature T c : The regions I, II, III in Fig. [2] are 
characterized by e = ln(T/T c ) <S 1 and h = H/H c2 (0) <sc 1. 
These are domains in which only classical AL fluctuations are 
important and the DOS contribution is less singular. The AL 
contribution to f} xtJ reads 



fS* = 2/3 Q F(x)/x, 



(11) 



where x = e/2h, fio = kse/nh = 6.68 nA/K is the quantum 
of thermoelectric conductance, and the function F(x) is given 
by Eq. (|3). One can see that the functional dependence of 
the corresponding contribution to the Nernst coefficient on x 
coincides with the qualitative estimate (0. The magnetization 
contribution to the observable /3 xy [see Eq. (0]i] is given by 



= A)[ln 



r(i/2 + jc) 

V2^ 



-xi//(l/2 + x) + x\. (12) 
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quantum 

IV quantum-to 
y y classical 
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FIG. 2: (Color online) Different asymptotic regions for the fluctua- 
tion Nernst coefficient on the H - T phase diagram. 



In the limit of vanishingly small magnetic fields h <K e (re- 
gion I), we find p xy = Bo(h/2e), which is two times larger 
than the result of Refs. flfl m. 221. The origin of the addi- 
tional factor should be traced back to the complicated analytic 
structure of the heat-current block (0. Thus the magnetiza- 
tion contribution/?^ = -fio(h/6e) cancels only 1/3 of fi xy , and 
the final result appears to be four times larger than the result 
of Refs. 03 03 El: 



neDH 
l2c(T-T c y 



k<e«l. 



(13) 



In the limit e <K h (region II), and close to the transition 
line, at h + e <sz h (region III), Eqs. (fTTT i and (fT2l yield 



^= J 8 [l-(ln2)/2] ! 



e + h 



H c2 (T) 
H-H c2 (T)' 



e«/i«l; (14) 
e + h « h « 1. (15) 



Now we turn to the quantum regime close to the upper crit- 
ical field H c2 (0) = ncT c /2yeD (regions IV, V, VI in Fig.©, 
where y = 1 .78 . . . is the Euler constant. Here the role of 
magnetization term becomes crucial: The 1/T divergence of 
= cM z /T exactly cancels the divergent contribution origi- 
nating from p xlJ , which is necessary to satisfy the third law of 
thermodynamics. As a result, the total coefficient f3 X y remains 
finite in the limit or zero temperature. The exact analytical ex- 
pression at t = T/T c «: 1 and rj = (H - H c2 (t)) / H c2 (t) «: 1 is 
quite lengthy. We present below only the asymptotic expres- 
sions in the regions IV, V, VI. 

In the purely quantum limit of vanishing temperature and 
away from H c2 (Q) (t <sz rj, region YV),fi xy is negative: 



UJ 2/3 () yt 



PqticT 



9eD[H - H c2 (0)] ' 



(16) 



This change of sign in thermoelectric response is similar to 
negative magnetoresistance in the quantum fluctuation trans- 
port for the usual electrical conductivity found in Ref. 12711 in 
the vicinity of H c2 (Q). The sign change is due to the DOS con- 
tribution being numerically larger than the positive AL term. 

In the quantum-to-classical crossover region, where H 
tends to H c2 {t) but remains limited as f 2 /ln(l/f) «; 77 «C t 
(region V), the coefficient f} XIJ becomes positive: 



/3 x Y y = £0 ln(r/?7), r / ln( 1 /r) « 77 « t « 1 . 
Near // c2 (f) (?7 « f 2 /ln(l/r), region VI), we find: 



?or 2 f 2 /3/7, 



77 « f 2 /ln(l/f) « 1. 



(17) 



(18) 



We also address the full classical region just above the tran- 
sition line, which covers a wide range of temperatures and 
magnetic fields (77 «: 1, region VII). In this region, f} xy is gen- 
erally comparable to = — /3o / rj, and we obtain 



PVTI - _ 



1 + 



h f"(l/2 + /i/4yf) 



4yt t//'(l/2 + h/4yt) 



rj^O, (19) 
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FIG. 3: (Color online) Comparison with experiment. Circles: ex- 
perimental data for lim H ^ f3 xy / H vs. e = lnT/T c obtained for the 
12.5-nm-thick Nbo.15Sio.85 film fioll . Dashed line: theoretical predic- 
tion for the strictly 2D geometry. Solid line: theoretical prediction 
for the real sample fioll with the 2D-3D crossover taken into account. 



with h = H/H c2 {0). Close to T c , Eq. (O matches Eq. ( fl5l . 
while in the limit T — > it matches Eq. ( fT8l ) provided that 
77 <§: t 2 1 ln(l/f). We note that in deriving Eq. 1191 . the account 
for Landau quantization of Cooper pair motion was crucial. 

Finally, we address the "non-singular" region far from the 
transition line (regions VIII and IX in Fig. [2j. Similar to the 
fluctuation correction to conductivity, the Kubo contribution 
j3 xy diverges as [lnln(l/T c T) - lnlnmax(/i, ?)], with 1/r play- 
ing role of the ultra-violet cutoff of the cooperon modes. Re- 
markably, exactly the same divergence of the opposite sign is 
contained in the magnetization contribution p^. Therefore, 
the total expression for /3 xlJ in this region is convergent and 
decreases nearly as a power law when increasing T and H: 



£vm=A>; 



eDH 



vm H 6ncT\n(T/T c y 

xy nc T 

^ IX ~ P° 12eDHln[H/H c2 (0)Y 



(l,/z)«f; (20) 
(l,f) (21) 



We see that even far from the transition the fluctuation Nernst 
coefficient can be comparable or parametrically larger than the 
Fermi liquid terms. In fact, it is conceivable that in some ma- 
terials the Cooper channel contribution to thermal transport at 
low temperatures can be dominant even in the absence of any 
superconducting transition at all (e.g., if superconductivity is 
"hidden" by another order). 

In Fig. [3] we compare our theory with the experimental data 
iflfil on the Nernst coefficient for the Nbo .isSio.ss film of thick- 
ness d = 12.5 nm. The dashed line shows our prediction for 
lim#^o {Fy/Hin a wide range of temperatures up to 30 T c (in- 
cluding regions I and VIII). We used the diffusion coefficient 
D = 0.087 cm 2 /s which is 60% of that reported in Ref. H3l 
(with k F l ~ 1, the precise determination of D is questionable). 
The data can be described by the 2D theory only in a close 
vicinity of T c , at e < 2. For larger e, the superconducting co- 
herence length £(T) becomes shorter than d and 3D nature of 
diffusion in the film should be taken into account. This can be 



done by substituting a„ — » a„ + D{npld) 2 and performing an 
additional summation over p = 0, 1, ... in Eqs. (HJ-O). The 
resulting curve is shown in Fig.[3]by the solid line. 

In summary, we have developed a complete microscopic 
theory of the fluctuation Nernst effect in a two-dimensional 
superconductor. Our results provide a natural explanation for 
a large Nernst signal observed in superconducting films Jiol 
1 111 and probably should be relevant to cuprates, where the 
energy scale of the chemical potential of preformed Cooper 
pairs is set by the pseudogap temperature T*. Our theoreti- 
cal predictions include a near-power decay of the transverse 
thermoelectric response away from the transition line, which 
is expected to persist well into the metallic phase. 
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